Decoherence induced by interacting quantum spin baths 
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We study decoherence induced on a two-level system coupled to a one-dimensional quantum spin 
chain. We consider the cases where the dynamics of the chain is determined by the Ising, XY , or 
Heisenberg exchange Hamiltonian. This model of quantum baths can be of fundamental importance 
for the understanding of decoherence in open quantum systems, since it can be experimentally 
engineered by using atoms in optical lattices. As an example, here we show how to implement a 
pure dephasing model for a qubit system coupled to an interacting spin bath. We provide results 
that go beyond the case of a central spin coupled uniformly to all the spins of the bath, in particular 
showing what happens when the bath enters different phases, or becomes critical; we also study 
the dependence of the coherence loss on the number of bath spins to which the system is coupled 
and we describe a coupling-independent regime in which decoherence exhibits universal features, 
irrespective of the system-environment coupling strength. Finally, we establish a relation between 
decoherence and entanglement inside the bath. For the Ising and the XY models we are able to 
give an exact expression for the decay of coherences, while for the Heisenberg bath we resort to the 
numerical time-dependent Density Matrix Renormalization Group. 
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I. INTRODUCTION 

Decoherence refers to the process through which super- 
positions of quantum states are irreversibly transformed 
into statistical mixtures. It is due to the unavoidable 
coupling of a quantum system with its surrounding envi- 
ronment which, as a consequence, leads to entanglement 
between the system and the bath. This loss of coherence 
may considerably reduce the efficiency of quantum infor- 
mation protocols [Ij] and it is crucial in describing the 
emergence of classicality in quantum systems 

Although desirable, it is not always possible to fully 
characterize the bath. Therefore it is necessary to resort 
to ingenious, though realistic, modelizations. Paradig- 
matic models represent the environment as a many-body 
system, such as a set of bosonic harmonic oscillators [3, |3| 
or of spin- 1/2 particles 0]. In some cases it is also pos- 
sible to recover the effects of a many-body environment 
via the coupling to a single-particle bath, provided its dy- 
namics is chaotic Q. In order to grasp all the subtleties 
of the entanglement between the system and its envi- 
ronment, it would be of great importance to study engi- 
neered baths (and system-bath interactions) that can be 
realized experimentally. In this paper we discuss a class 
of interacting spin baths which satisfy these requirements: 
a two-level system coupled to a one-dimensional array of 
spin-1/2 particles, whose free evolution is driven by a 
Hamiltonian which embraces Ising, XY and Heisenberg 
universality classes. In several non-trivial cases we can 
solve the problem exactly. Moreover we show that it is 
possible to realize these system-l-bath Hamiltonians with 
cold bosons in optical lattices. We thus extend the ap- 
proach of Jane et al. Q showing that optical lattices can 
be useful as open quantum system simulators. 

Our analysis can be framed in the context of the re- 



cently growing interest in the study of decoherence due 
to smn baths, see [1, i, 0, [U, [H, [H, [11, [H, IH, 
[isl . [igl . [20j . Starting from the seminal paper of Zurek 
several papers analyzed the decoherence due to a collec- 
tion of independent spins. Cucchietti, Paz and Zurek Q 
derived several properties of the Loschmidt echo start- 
ing from fairly general assumptions for the distribution 
of the splittings of the bath spins. The effect of an 
infinite- rang e interaction among the spins was introduced 
m Rcf. [lOl; the same model was further exploited in 
Dawson et al. [ll| to relate decoherence to the entangle- 
ment present in the bath. Properties of decoherence in 
presence of symmetry breaking, and the effect of critical 
behavior of the bath were discussed in [l^, Most 
of the works done so far are based on the so called cen- 
tral spin model, where the two-level system is coupled 
isotropically to all the spins of the bath. This assump- 
tion tremendously simplifies the derivation, but at the 
same time it may introduce some fictitious symmetries 
which are absent in realistic systems. Moreover it can be 
very hard to simulate it with engineered baths. A crucial 
feature of our work is the assumption that the two-level 
system interacts with only Jew spins of the bath. As we 
will show, this introduces qualitative differences as com- 
pared to the central spin model, moreover it is amenable 
to an experimental implementation with optical lattices. 

The paper is organized as follows. In the next Section 
we introduce the system-|-bath model Hamiltonians that 
will be studied throughout this paper and we define the 
central quantity that will be used in order to character- 
ize the decoherence of the system, i.e. the Loschmidt 
echo. We then show in Section IIIII how optical lattices 
can be used to simulate the class of Hamiltonians intro- 
duced in Section |TT1 The rest of the paper is devoted to 
the derivation and to the analysis of our results, both 
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for a single system-bath link (Section lIVp . and for mul- 
tiple links (Section |V]). When the two-level system is 
coupled to an XY- or XX-model, it is possible to de- 
rive an exact result for the Loschmidt echo. This is ex- 
plained in Section flV Al where we also discuss in detail 
its short- and long-time behavior and relate it to the crit- 
ical properties of the chain. Further insight is obtained 
by perturbative calculations which agree very well, in 
the appropriate limits, with the exact results. In Sec- 
tion |TVB] we present our results for the Heisenberg bath. 
In this case an analytic approach is not possible. Here 
we solve the problem by means of the time-dependent 
Density Matrix Renormalization Group (t-DMRG - see 
Appendix |A| . In Section IIV CI we analyze the possible 
relation between decoherence and entanglement proper- 
ties of the environment: we relate the short-time decay of 
coherences to the two-site nearest-neighbor concurrence 
inside the bath. In Section |V] we extend our results to 
the case in which the system is coupled to an arbitrary 
number of bath spins; a regime in which the decoherence 
is substantially independent of the coupling strength be- 
tween the system and the environment is discussed in 
Section IV Al Finally, in Section IVII we draw our con- 
clusions. In the Appendices we give some technical de- 
tails on the numerical DMRG approach (Appendix 
we provide explicit expressions for the Fermion correla- 
tion functions needed for the evaluation of the Loschmidt 
echo in the XF- model (Appendix [B]) and we briefly dis- 
cuss the central spin model (Appendix [C| to make our 
paper self-contained. A brief account of some of the re- 
sults discussed in this paper is presented in [2l| . 



II. THE MODEL 

The model we consider consists of a two-level quantum 
object (qubit) S coupled to an interacting spin bath E 
composed by N spin-1/2 particles (see FigH^): the idea 
is to study how the internal dynamics of E affects the 
decoherent evolution of S. 

In our scheme the global system -t- i? is fully charac- 
terized by a standard Hamiltonian of the form 

H = Hs + He + 'Hint , (1) 

with Hs,E being the free-Hamiltonians of S and E, and 
Hint being the coupling term. Without loss of generality 
we will assume the free Hamiltonian of the qubit to be 
of the form 

Ws = y (l-r.)=c.,|e)(e| , (2) 

with Tq, being the Pauli matrices of 5 (a = x,y,z), and 
|e) its excited state (the ground state being represented 
by the vector \g)). On the other hand, the environment 
will be modeled by a one-dimensional quantum spin-1/2 



chain described by the Hamiltonian 

i 

+ Aa|a|+i + 2Aa|] , (3) 

where af (a = y, z) are the Pauli matrices of the 
i-ih spin. The sum over j goes from 1 to — 1 for 
open boundary conditions, or from 1 to for periodic 
boundary conditions (where we assume that (J%j^i = a"). 
The constants J, A, 7 and A respectively characterize 
the interaction strength between neighboring spins, the 
anisotropy parameter along z and in the xy plane, and an 
external transverse magnetic field. The Hamiltonian Q 
has a very rich structure [22| . For the sake of simplicity, 
we shall consider the following paradigmatic cases: 

• the XY -model in a transverse field - see Sec. IIV Al 
Here one has A = and A, 7 generic. For < 7 < 
1, Eq. ([U belongs to the Ising universality class, 
and it has a critical point at |Ac| = 1; for 7 = 
it reduces to the XX universality class, which is 
critical for |A| < 1. 

• the XXZ anisotropic Heisenberg model - see 
Sec. IIVBI Here one has A,7 = and A generic. 
In this case the Hamiltonian ^ is critical for 
— 1 < A < 1 while it has ferromagnetic or anti- 
ferromagnetic order for A>lorA<— 1 respec- 
tively. 

Finally, the qubit S is coupled to the spin bath through 
a dephasing interaction of the form 

Wint = -e £ |e) (ek| , (4) 
j=h 

where e is the coupling constant, and the link number 
m counts the number of environmental spins (labeled by 
ji ■ ■ -jm) to which 5" is coupled — FiglT^ refers to the 
case where S is interacting with the first spin of an open- 
boundary chain, i.e. m — 1, ji = 1. 

By varying the parameters m, A, 7, A and e, the above 
Hamiltonians allow us to analyze several non-trivial S+E 
scenarios. Moreover we will see in Sec. lIIII that it is possi- 
ble to use optical lattices manipulation techniques [3, [2^ 
to experimentally simulate the resulting dynamical evo- 
lution. 



A. Decoherence: the Loschmidt echo 

With the choice of the coupling ^ the populations of 
the ground and excited states of the qubit do not evolve 
in time, since [t^,H] = 0. Consequently, no dissipation 
takes place in the model and the qubit evolution is purely 
decoherent: the system S looses its coherence without 
exchanging energy with the bath [U, [2^ . To study this 
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effect of pure phase decoherence, we suppose that at the 
beginning the qubit is completely disentangled from the 
environment — namely, at time t — the global system 
wave-function is given by 

\^{Q))^\m)s'»\m)E, (5) 

where |(/)(0))s = Cg \g)+Ce \e) is a generic superposition of 
the ground and the excited state of S and |</3(0)) e = \G)^ 
is the ground state of the environment Hamiltonian TIe- 
The global evolution of the composite system under the 
Hamiltonian Eq. ([T]) will then split into two terms: one 
where E evolves with the unperturbed Hamiltonian Tig = 
Ti.E and the other where E evolves with the perturbed 
Hamiltonian 

He = He + {e\nint\e) . (6) 

As a result, one has 

|vI/(0)) ^ |vl'(t)) = cg\g)^\^g{t))^ (7) 

+ Cee-*"=*|e)(^|^e(i))B 

where the two branches of the environment 
are Mt)) = e-*«»* 1^(0))^ and \Mt)) e = 
g-'in^t 1^(0))^. Therefore the evolution of the re- 
duced density matrix p = Tie {'^\ of the two-level 
system corresponds to a pure dephasing process. In the 
basis of the eigenstates {\g), |e)}, the diagonal terms pgg 
and Pee do not evolve in time. Instead the off-diagonal 
terms will decay according to 

Peg{t)^Peg{0)e-'^^'D{t) , 

where 

Dit) EE {^git)\^eit)) = (^(0)|e'««*e-'«=*|¥>(0)) (8) 

is the decoherence factor. The decoherence of S can then 
be characterized by the so-called "Loschmidt echo" , i.e. 
by the real quantity 

C{t) = = |(G|e-*(«^+<^l«'"'l^»*|G)|2 , (9) 

where in the equality term we used the fact that |'/'(0)) = 
\G) is the ground state of He- On one hand, values 
of £{t) close to 1 indicate a weak interaction between 
the environment and the qubit (the case C{t) = 1 cor- 
responds to total absence of interaction). On the other 
hand, values of £{t) close to correspond instead to a 
strong suppression of the qubit coherence due to the in- 
teraction with E (for C{t) = the qubit is maximally 
entangled with the environment, and its density matrix 
p becomes diagonal). In the following we will analyze 
the time evolution of the Loschmidt echo to study how 
the internal dynamics of the environment E affects the 
decoherence of the system S. We will start with the case 
when S is coupled to a single spin of the chain E (i.e. 
m = 1, see Sec. IIV[) . The case with multiple links will be 
instead discussed in Sec. |Vl 



For the sake of completeness, we finally notice that, 
since the function © measures the overlap between the 
time-evolved states of the same initial configuration un- 
der two slightly different Hamiltonians, one can use it as 
an indicator of the stability of motion. This kind of anal- 
ysis has been performed in Ref. 26], where the fidelity 
between the ground states of quadratic Fermi Hamiltoni- 
ans has been analyzed in connection with quantum phase 
transitions. 

System EnYironment 

TzUi. 




FIG. 1: a) A sketch of the system-plus-bath model we con- 
sider in this work. The two-level system (at position zero) is 
coupled to the a" component of the first spin of the chain that 
acts as a spin bath. Atoms in an optical lattice can simulate 
this controlled decoherence by means of series of lasers b) and 
lattice displacements c), which allow to realize both the in- 
teraction of the bath with an external magnetic field and the 
anisotropic exchange coupling present in Eq. ((3}. 



III. SIMULATION OF OPEN QUANTUM 
SYSTEMS BY OPTICAL LATTICES 

Before analyzing in detail the temporal evolution of the 
Loschmidt echo (0) , in this Section we present a method 
which would allow one to experimentally simulate the 
dynamics induced by the Hamiltonian of Eq. ([1]) in a 
realistic setup. 

The S + E system introduced in Sec. [TTl can be seen 
as an "inhomogeneous" spin network with -I- 1 sites, 
where one of the spins (say, the first) plays the role of the 
system of interest S while the remaining N play the role 
of the environment E. This immediately suggests the 
possibility of simulating the dynamical evolution of such 
a system on optical lattices by employing the techniques 
recently developed in Refs. 0, HI] . An important aspect 
of our scheme however is the fact that we assume the 
coupling between the system S and E to be independent 
from the couplings among the N spins which compose the 
environment. Analogously the free Hamiltonian of S is 
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assumed to be different from the on-site terms of the free- 
Hamiltonian of E. On one hand this aUows us to study 
different environment Hamiltonians without affecting the 
coupling between E and S. On the other hand, this also 
allow us to analyze different S + E coupling regimes (e.g. 
strong, weak) without changing the internal dynamics of 
the bath. 

In order to include these crucial ingredients we found 
more convenient to follow the scheme developed by Jane 
et al. in Ref. (the optical lattice manipulation schemes 
developed by Duan et al. in Ref. [1^ seem to be less ad- 
equate to our purposes). The key advantage is that we 
can realize the system-plus-bath setup by using a single 
one-dimensional lattice in which the quantum system is 
placed on a given lattice site (for example the first one, 
as in Fig.[T}D). The different Hamiltonians for the system 
and for the bath are realized by specific pulse sequences 
which simulate the dynamics of the model. The same 
holds for the coupling Hamiltonian of the two-level sys- 
tem with the bath, which is different from the couplings 
within the bath. In Fig. [1] the leftmost atom simulates 
the two-level system, the coupling to the second site is 
the interaction between the quantum system and the en- 
vironment, the rest of the chain is the interacting spin 
environment. 

Jane et al. 0] showed that atoms loaded in an opti- 
cal lattice can simulate the evolution of a generic spin 
Hamiltonian in a stroboscopic way when subjected to 
appropriate laser pulses, Fig.IIh, and controlled displace- 
ments. Fig. [TJ;, which allow to implement the single-site 
and two-site contributions to the Hamiltonian. The key 
point is that in our case the sequences of gates need to al- 
low for discriminating between the system and the bath. 
The types of baths that one can simulate by these means 
embrace Ising, XY and Heisenberg exchange Hamilto- 
nian. Therefore, by varying the parameters of the op- 
tical lattice, it is possible to test the impact of the dif- 
ferent phases (critical, ferromagnetic, anti-ferromagnetic, 
etc. . . ) of the environment on the decoherence of the two- 
level system. 

Following the idea of a Universal Quantum Simulator 
described in Q, the time evolution operator associated 
with 7i over a time t can be simulated by decomposing 
it into a product of operators acting on very short times 



e-*^* = hm 



N 



(10) 



TJ-XX / J(1 + 7)t '| J-J-yy ( .7(1 'I)t \ jjzz ( J^T \ 

^Jj+U 2 y^ij + U 2 j '^Jj+U 2 / 



where r = t/n, U^{e) = e*^'^^ Ul^{0) = e'^"'^"'-, 
and the index labels the two-level system S. For 
a e {x,y} one can write Uf^ = V°'V^U'r^^y^^V°'\ 

where 1/" = (1 - icr")/\/2 are fast homogeneous local 
unitary operations. These can be realized with single 
atoms trapped in an optical lattice 0], each having two 



relevant electronic levels (|0)j , |l)j) interacting with a 
resonant laser according to: 



= n{e"<'\i) 



|o),(i|) 



(11) 



The evolution under the Hamiltonian Eq. (|lip . 



U^it, 



, yields the single-qubit operations 



TLrjT 



UH4r,0)UHMA) and 



i ^20' 



f^f(^) = '^/(^''^+^)t^, (#T:0): whence C/,n-f) 



while Uj"{^) — i(Jj. These operation can be made very 
fast by simply increasing the laser intensity and thereby 
the Rabi frequency. They can be performed either si- 
multaneously on all qubits, by shining the laser homoge- 
neously onto all atoms, or selectively on some of them, by 
focusing it appropriately (see Fig.[T|D). For our purposes 
the individual addressing is needed only for the atom in 
position 0, which represents the quantum system; this is 
anyway the minimal physical requirement for being able 
to monitor its state during the evolution. 

Two-qubit operations can be performed by displacing 
the lattice in a state-selective way [23|, so that state 
|0)j|l)j+i acquires a phase factor e~"^, as experimen- 
tally realized in [2^. The resulting gate Gj,j+i{ip) can 
be composed with rotations to yield Ujj_^_i{6) — 
e''^[Gjj+i{2e) cr^a"^^^]^. This will affect all atoms from 
to N. Since we want a different coupling for the {01} pair 
than for all others, we need to erase the effect of the inter- 
action for that specific pair using only local operations on 
atom 0, as in the sequences [cg U§i{9)]'^ — [ctq J/gf (6*)]^ = 
K USnO)? = a. Defining C/|"(0) ^ JJU C^""+i(^), we 
can generate each simulation step in Eq. (|10p as 



^ [a-C/-((e-^)i)]'K-L^-(- 



J(7+1)tN]2 



)] (12) 



involving only global lattice displacements, global laser- 
induced rotations and local addressing of atom 0. 

We note that an alternative scheme exists [l^l, based 
on tunnel coupling between neighboring atoms rather 
than on lattice displacements, which can attain the simu- 
lation described here for the special case 7 = and would 
require some additional stroboscopic steps in order to re- 
produce the general case. 

Apart from ID spin baths, this approach can be ex- 
tended to other types of environment. For example, it 
would be quite interesting to consider, as an engineered 
bath, a 3D optical lattice. Besides being feasible from 
an experimental point of view, this could be useful in 
studying for instance the situation found in solid-state 
NMR . It would be also intriguing to study the Bose- 
Hubbard model as a bath, which would make the exper- 
imental realization even simpler [sot . Here we just focus 
on one-dimensional baths since, in several cases, they are 
amenable to an exact solution. 
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IV. RESULTS: THE SINGLE-LINK SCENARIO 



where 



: for Tia, while 



In this Section we analyze the time evolution of the 
Loschmidt echo from Eq. Q for several distinct scenar- 
ios where the qubit S is coupled to just one spin of the 
chain, i.e. m = 1 in Eq. (|4]). In this case, in the thermo- 
dynamic limit the interaction Tiint between the system 
S and the environment E does not affect the description 
of the bath Hamiltonian He, since it is local; therefore 
it can be considered in all senses as a small perturbation 
of the environment. The bath is effectively treated as a 
reservoir, which is in contact with the system through 
just one point. In the following we study the cases in 
which the bath is described by a one dimensional spin- 
1/2 Ising, an XY (Subsec. lIV Ap and a Heisenberg (Sub- 
sec. |IVB]) chain. 



A. XY-hath 

Here we focus on the case of a spin bath E charac- 
terized by a free Hamiltonian ^ of the XY form, i.e., 
with null anisotropy parameter along z (A = 0). In 
this case the environment is a one-dimensional spin-1/2 
XY-mode\, which is analytically solvable [31]. Below we 
show that also the Loschmidt echo can be evaluated ex- 
actly [2l| . both for open and for periodic boundary con- 
ditions. 

In the case m < in which the system S is cou- 
pled to just some of the spins of the bath E, the per- 
turbed Hamiltonian He of Eq. © is the Hamiltonian of 
an XY chain in a non-uniform magnetic field. In this cir- 
cumstance one cannot employ the approach of Ref . , 
and in general the dynamical evolution of the system 
has to be solved numerically. The derivation we present 
here instead is analytical and it applies for all values of 
m — 1, ■ ■ ■ , N. In the following we will present it for the 
case of a generic m but, in the remaining of the Section, 
we will explicitly discuss its results only for the single- link 
case (i.e. m = 1). 

The first step of the analytical derivation is a Jordan- 
Wigner transformation (JWT), in order to map both 
Hamiltonians Hg and He onto a free-Fermion model, de- 
scribed by the quadratic form [3l| 



1 



9 X! "^^'^ ' 



(13) 

where Ci,c| are the annihilation and creation operators 
for the spinless Jordan- Wigner fermions, defined by 

Cfe = exp (^in ^ cr+(T^"^ . 

The two matrices A, B are given by 

= -J(4.j+i+^j,fc+i)-2(A + e,)<5,,fe, (14) 



(15) 



if S is coupled to the j-th spin 
elsewhere 



(16) 



for He- A generic quadratic form, like Eq. p3|) (where 
A is a Hcrmitian matrix, due to the Hermiticity of H, 
and B is antisymmetric, due to the anti-commutation 
rules among the Ci), can be diagonalized in terms of the 
normal- mode spinless Fermi operators {ilk,Vk}' 



H = Y,Ek 



1 



(17) 



where rik = J2ii9k,tCi + hk,ic\), or in matrix form: 

j7=g.c + h-c^. (18) 

If we rewrite the two change-of-basis matrices g and h as 
9kA = \{<i>k,i + i'k,i) and hk,i = \{4ik,i - ipk,i), we even- 
tually arrive at the following coupled linear equations, 
whose solution permits to find the eigenbasis of the non- 
uniform Hamiltonian in Eq. (I13p : 



0fc(A- 



B) = 
B) = 



Eki^k, 
Ek4>k- 



(19) 



Since A is symmetric and B is antisymmetric, all of the 
Ek are real; also the gk,i and the hk,i can be chosen to be 
real. The canonical commutation rules for the normal- 
mode operators impose the constraints: g g"^ + h h"^ = 1; 

gh^-hgT = 0. 

It is convenient to rewrite the spin-bath Hamiltonian- 
plus-interaction He in the form 



He = i*^C* 



(20) 



where — . . . cjy Ci . . . c^^ [ci are the correspond- 
ing spinless Jordan- Wigner Fermion operators) and C = 
(T^ (g) A -I- ia^ (g) B is a tridiagonal block matrix. 

The Loschmidt echo, Eq. can then be evaluated 
by means of the following formula [13] : 



= |det(l-r + re-'C*)| , 
(21) 

where the elements of the matrix r are simply the two- 
point correlation functions of the spin chain: rij — 

{G\^\^j\G), where \G) is the ground state of Hg. An 
explicit expression of the correlators as well as of the 
matrix e~**^* is given in Appendix [Bl 



Equation ([2T|) provides an explicit formula for the 
Loschmidt echo in terms of the determinant of a 2N x 2A^ 
matrix, whose entries are completely determined by 
the diagonalization of the two linear systems given by 
Eq. . This is one of the central results of our work. 
It allows us to go beyond the central spin model where 
all the spins of E are uniformly coupled with S (m = N)^ 
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whose solution (at least for periodic boundary condi- 
tions) was discussed in [l^ [fl [IB| and, for the sake 
of completeness, has been reviewed in Appendix \C\ (see 
Eq. (jC6[) ). Notice that, similarly to the central spin 
model, this formula allows to study a system composed 
of a large number of spins in the bath N ^ 10^ — 10"^, 
as it only requires manipulations of matrices whose size 
scales linearly, and not exponentially, with N. 

1. Ising bath: general features 

The generic behavior of £ as a function of time for 
different values of A, and fixed coupling constant e, is 



shown in Fig. [2] For A < 1 the echo oscillates with a fre- 
quency proportional to e, while for A > 1 the oscillation 
amplitudes are drastically reduced. The Loschmidt echo 
reaches its minimum value at the critical point Ac = 1, 
thus revealing that the decoherence is enhanced by the 
criticality of the environment. Since the chain is finite, 
at long times there are revivals of coherence [SS^; in the 
thermodynamic limit iV — > oo they completely disappear. 
In any case, as it can be seen from the figure, already for 
N = 300 spins there is a wide interval where the asymp- 
totic behavior at long times can be analyzed. 




FIG. 2: Single-link Ising model (m = 1, A = 0,7 = 1): Loschmidt echo of Eq. (I21|l as a function of time for a qubit coupled 
to a = 300 spin chain with periodic boundary conditions (here e = 0.25). The various curves are for different values of the 
transverse magnetic field: A = 0.25 (black dotted), 0.5 (red), 0.9 (green), 0.99 (blue), 1 (brown), 1.01 (violet), 1.1 (magenta), 
1.5 (orange), 1.75 (cyan). The critical point corresponds to A = 1 (brown curve). Notice the revivals of quantum coherence at 
t* ~ 150, due to the finite size of the chain [3^ . 



2. Short-time behavior 

At short times the Loschmidt echo C decays as a Gaus- 
sian [lil: 

Cit) - e""*' , (22) 



as it can be seen in Fig. [3^, which shows a magnification 
of the curves from Fig. [5] for small times. This behavior 
can be predicted within a second-order time perturbation 
theory in the coupling e between the system and the bath: 
if e is small as compared to the interaction J between 
neighboring spins in the bath (e <C J), the decoherence 
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0.994 
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FIG. 3: Short-time behavior of the Loschmidt echo: a) mag- 
nification of the plot in Fig. [2] for t < 0.6, where a Gaussian 
decay L ~ e~°* is visible; b) dependence of a on A (the solid 
red line shows the perturbative estimate given by Eq. (|26|) V 



factor D{t) of Eq. ^ can be expanded in series of e: 

ft 



T 



exp I ~i I dt'e^'^^^'We-''^^^' 



1 + eAi +6^X2, 



(23) 



where T is the time ordered product and W = e erf ac- 
counts for the interaction of the two-level system with 
the spin chain. The above expression has to be evalu- 
ated on the ground state of Hg, therefore it is useful to 
rewrite the interaction W in terms of the normal mode 
operators r][^^ of Hg: 



The first-order term then reads 



(24) 



while the second-order term is given by 



/ dt' f dt" 4V((g,,i%i) 
"'o Jo L , /„■ ^ 



i(E,+Ej){t'-t") 



(25) 



The Loschmidt echo is then evaluated by taking the 
square modulus of the decoherence factor: 

C{t) ~ 1 - 462^2 ^ [(5,,i/i,,i)2 - 5,43,- 1] . (26) 




FIG. 4: Single-link Ising model (m = 1, A = 0,7 = 1): Be- 
havior of the Loschmidt echo of Eq. (|21|l at short times. The 
plot shows the rescaled parameter dxoi/e^ as a function of A 
for a periodic A'^ = 200 Ising spin chain. Various symbols are 
for different values of the coupling strength e. Inset: plot in 
a semi-logarithmic scale. The dashed line indicates a fit of 
numerical data; d\oi/(^ = ci log |AA| -f C2, with ci ~ 0.40983 
and C2 ^ 0.000108. 



In Fig. the initial Gaussian rate a is plotted as a 
function of A; circles represent numerical data, while the 
solid red curve is the perturbative estimate obtained from 
second-order perturbation theory, given by Eq. (j26p . In 
Fig. m we analyze the behavior of the first derivative of the 
rate a{\, e) as a function of the distance from criticality 
AA = A — Ac, for a fixed number N of spins in the chain. 
As predicted by the perturbative estimate, a scales like 
e^; most remarkably, its first derivative with respect to 
the transverse field diverges if the environment is at the 
critical point Ac- In the inset we show that d\a diverges 
logarithmically on approaching the critical value, as: 



da 



ci In |A — Ac| + const. 



(27) 



This is a universal feature, entirely due to the underlying 
criticality of the Ising model. 

Our results show that at short times the Loschmidt 
echo is regular even in the presence of a bath undergoing 
a phase transition. The critical properties of the bath 
manifest in the changes of C when the bath approaches 
the critical point. 



3. Long-time behavior 

At long times, for A > 1 the Loschmidt echo ap- 
proaches an asymptotic value £00, while for A < 1 it 
oscillates around a constant value (see Fig. [5] for a qual- 
itative picture). This limit value £00 strongly depends 
on A and presents a cusp at the critical point, as shown 
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FIG. 5: Single-link Ising model (m = 1, A = 0, 7 = 1): Long- 
time behavior of the Loschmidt echo. The two plots show the 
saturation value Lac as a function of A in an A*' = 200 Ising 
spin chain with periodic b.c, for two different values of the 
coupling strength: e — 0.05, 0.25. Circles indicate numerical 
data, solid line is obtained from the perturbative calculation 
in Eq. 



in Fig. [5l Evidence that £00 describes the asymptotic 
regime can be obtained by comparing data with the re- 
sult of an analytical expression based on the following 
simple ansatz: We assume that at the thermodynamic 
limit N 00 the coherence loss saturates and is con- 
stant after a given transient time to 



yt>t 



0- 



(28) 



By writing the ground state \G) of Tig on the eigenbasis 
of Tie, the Loschmidt echo in Eq. ^ can be formally 
written as 



C 



(29) 



Here, Ck are the coefficients of the expansion |'0o^'') = 



E,c.|^i^)), where He\4'^) = 
|G). An integration of the ansatz Eq 
T > to to time T + A thus gives: 



E 



Cfe e 



. 

/ smc — 



and l^'')) ^ 
from time 



(30) 



where sine (x) — sin(a;)/a;. If the ground state in the per- 
turbed Hamiltonian Tie is not degenerate, we can just 
retain the term in the sum corresponding to A: = 0, 
thus obtaining Coo ~ |co|*- Therefore a rough estimate 
of the Loschmidt echo is given by the overlap between 
the ground state of Hg and the ground state of He (see 



Ref. ^2d\ for a detailed analysis of the ground state fidelity 
in spin systems) . This can be evaluated via second-order 
perturbation theory in e: 



Co 



-E 



Ei'^ - E. 



(a) 



(31) 



where 



E 



(a) 



are the excited states of Tig with energies 

and W — e X]j=ji ■ The expression in Eq. ([3T|) 
can be easily computed for the case m = 1 following the 
same steps used to evaluate Eq. The final result is: 



l-2e' 



E 



gn h 



E, 



ia) 



E. 



(a) 



-\ 4 



(32) 



where the energies E^: 
energy, i?^^-* 



(9) 



_ j^ia) _ j^ia) 



are rescaled on the ground state 

In Fig. [5] we plotted these re- 
sults for two different values of the perturbation strength 
(solid lines); they can be compared with the exact nu- 
merical results (circles). Notice that analytic estimates 
are more accurate for small e, and far from the critical 
region. 



4-. Ising bath: scaling with the size 

As one could expect, contrary to the case of a non- 
local system-bath coupling , if the environment is not 
at the critical point, the decay of the coherences is inde- 
pendent of the number of spins in the bath. This is due 
to the fact that correlations away from criticality decay 
exponentially, while they decay as a power law at the 
critical point. In particular, we checked that both the 
initial rate a and the saturation value Coo do not depend 
significantly on N for A ^ Ac. We present here a simple 
argument as a confirmation of this statement. Far from 
the critical region, independently of the bath size (for 
sufficiently large baths), a qubit simultaneously coupled 
to two non-communicating iV-sites chains via a one-spin 
link (configuration 1) should decohere in the same way 
as if it was coupled to two sufficiently distant spins of a 
2A''-sites chain (configuration 2). A pictorial representa- 
tion of the two configurations described above is shown in 
Fig. [SI where the red spin represents the system, the grey 
ones are the bath spins directly coupled to the system. 

We numerically checked that, far from the critical re- 
gion, the two configurations leave the decoherence un- 
changed. In particular we chose an iV = 200 chain, and 
obtained that, if |A — Ad > 10"^, the Loschmidt echo 
evaluated in the two configurations is the same, up to a 
discrepancy of 10^^, as it can be seen from Fig. [T] 
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FIG. 6: Qubit (red) coupled to two non-communicating A'^- 
sites chains (configuration 1) and to two spins of a 2A''-sites 
cliain (configuration 2). Here, as an example, we set A'^ — 11. 
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FIG. 7: Upper figure: solid curves stand for the qubit coupled to two non communicating open A'' = 200 spin Ising chains 
(configuration 1); dashed curves are for the qubit coupled to two opposite spins of a periodic A'' = 400 spin Ising chain 
(configuration 2). Different curves are for various values of the transverse magnetic field: A = 0.5 (black), 0.9 (red), 0.99 
(green), 1 (blue), 1.01 (brown), 1.1 (magenta), 1.5 (cyan). Lower figure: absolute differences between the two configurations. 



At the critical point instead the behavior dramatically 
changes. In Fig. [5] we plot the decay of L in time at 
Ac for different sizes of the chain and fixed perturbation 
strength, e = 0.25. The minimum value £g reached by 
the Loschmidt echo decays with the size N as 



•^0 



l + /31n(iV) 



(33) 



Notice also that revivals are due to the finiteness of the 
environmental size (as in Fig. [5]) [sl] . We expect that the 
coherence loss should go to zero at the thermodynamic 
limit iV — > c»; this is hard to see numerically, since the 
decay is logarithmic, and the actual value of £q is still 
quite far from zero, even for 'N = 2000 spins. 



5. XY-baths with arbitrary xy-anisotropy 7 

The properties described in the previous Subsections 
are typical of the Ising universality class, indeed they 
remain qualitatively the same as far as A = 0,7 7^ 
in Eq. in particular this class of models has one 

critical point at A = 1. Following our previous analysis 
of the Ising model (7 = 1), we focus on the two distinct 
regions of short- and long-time behavior, whose features 
are depicted in Fig. [31 

At small times the decay is again Gaussian, as in 
Eq. (|^ . and the first derivative of the decay rate a 
with respect to A diverges at Ac (see Fig [5^). As far 
as the xy anisotropy decreases, the system approaches 
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FIG. 8: Single-link Ising model (m = 1,A = 0,7 = 1). 
Loschmidt echo as a function of time at the critical point 
A = 1, for different sizes of the periodic chain: A'' = 50 (black), 
100 (red), 200 (green), 300 (blue), 400 (brown), 500 (violet), 
1000 (orange), 2000 (cyan). The perturbation strength is kept 
fixed; e = 0.25. The dashed line is a guideline that shows a 
decay of the type ii'^(t) = co/(l + Cilnt). Inset: Minimum 
value of as a function of A'^. Numerical data have been 
fitted with LI = £o/(l + /31n7V) (dot-dashed line), where 
£0 ~ 0.99671, /3 ^ 2.36933 x 10"^ 
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FIG. 9: Single-link XF-models (m = 1,A = 0): Short- and 
long-time behavior of the Loschmidt echo for a qubit cou- 
pled to a periodic A'' = 300 spin XY chain, with a coupling 
strength e = 0.25. a) Small times (initial Gaussian decay 
rate q); b) Long times (saturation values £00). The various 
curves are for different anisotropy values: 7 = (XX-model, 
full black circles), 0.1 (red squares), 0.25 (green diamonds), 
0.5 (blue triangles up), 0.75 (brown triangles down), 0.9 (ma- 
genta triangles left), 1 (Ising, cyan triangles right). 



a limiting case in which d^a = 



a = —00. 



while d^c 



= 0. 



The asymptotic value Coo at long times shows a cusp at 
A = 1; notice that, as the anisotropy 7 decreases, de- 
spite the fact that the decay of coherences at short times 
is always reduced, their saturation value £00 becomes 



lower for A < 1, while gets higher for A > 1. Therefore 
decoherence is asymptotically suppressed only in the fer- 
romagnetic phase, while it is enhanced for A < 1. As we 
said before, when 7 — > the system is driven towards a 
discontinuity at A = 1. This is due to the fact that for 
A > 1 the ground state approaches a fully polarized fer- 
romagnetic state. Finally we analyzed the scaling of the 
decoherence with the size of the environment at critical- 
ity. As in the Ising model, the minimum value Cq reached 
by the Loschmidt echo at Ac depends on the chain length 
N as in Eq. = Co/{l + pin{N)). 




FIG. 10: Single-link XX-model (m = 1, A = 7 = 0): Tempo- 
ral behavior of the Loschmidt echo for a qubit coupled to a 
periodic A'^ = 300 spin isotropic XX chain; e — 0.25. The var- 
ious curves stand for different values of the magnetic field: A = 
0.25 (black), 0.5 (red), 0.75 (green), 0.9 (blue), 0.95 (brown), 
0.99 (magenta), 1. (orange). Bottom left: first derivative of 
the decay rate a at short times with respect to A, as a function 
of A. Bottom right: saturation value Lao at long times. 



The system behaviour dramatically changes when it 
becomes isotropic, i.e. at 7 = 0. In this case the envi- 
ronment is an XX spin chain and it does not belong any 
more to the Ising universality class. In particular it ex- 
hibits a critical behavior over the whole parameter range 
|A| < 1, while it is ferromagnetic (anti- ferromagnetic) 
for A > 1 (A < —1). The decay of the Loschmidt echo 
(shown in Fig. [TO)) reflects its critical properties: indeed 
we found that it behaves as in Eq. ([55)1 over the whole 
range |A| < 1. In the ferromagnetic case instead the 
coupled qubit does not decohere at all (i.e. C{t) = 1), 
since the ground state of the XX-mode\ for A > 1 is the 
fully polarized state with all spins parallel to the external 
field. The first derivative of the initial decay rate d\a is 
monotonically negative in the critical region by increas- 
ing A and diverges for A ^ 1~, while it is strictly zero for 
A > 1. Also the plateau £00 as a function of A presents a 
discontinuity in A = 1, since it drops to zero for A ^ 1~, 
and it equals one for A > 1. 
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B. XXZ-Heisenberg chain bath 

In this Section we consider a single-link (to — 1) XXZ 
anisotropic Heisenberg chain (A = 7 = in Eq. ^) with 
anisotropy parameter A. We resort to the numerical t- 
DMRG to compute the Loschmidt echo (O, since this 
specific spin model is not integrable. A brief introduction 
to the static and the time-dependent DMRG algorithm is 
given in Appendix [Xj We consider open boundary condi- 
tions, since DMRG with periodic boundaries intrinsically 
gives much less accurate results ,35]. The evaluation of 
the Loschmidt echo with the t-DRMG is more time- and 
memory-consuming than its exact calculation in the solv- 
able cases, therefore at present we could not study the 
behavior of at long times. Nonetheless we are able 
to fully analyze the short-time behavior in systems with 
environment size of up to ~ 10^ spins. In the fol- 
lowing we show the results concerning the coupling of 
the two-level system to one spin in the middle of the 
chain: this corresponds to the case with less border ef- 
fects. We numerically checked that the results are not 
qualitatively affected from this choice, but changing the 
system-environment link position results in a faster ap- 
pearance of finite size effects due to the open boundaries. 
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FIG. 11: Single-link XXZ-model (m = 1,A = 7 = 0): 
Loschmidt echo as a function of time for a qubit coupled to the 
central spin of an open ended A*' = 100 spin XXZ-Heisenberg 
chain, with coupling strength e = 0.1. The various curves are 
for different values of the anisotropy: A = 0.9 (black), 0.5 
(red), (green), -0.5 (blue), -0.9 (brown), -1 (magenta), - 
1.5 (violet), -2 (cyan), -2.5 (orange); curves corresponding to 
non-critical situations are dashed. 



In Fig.[TT]we plot the decay of the Loschmidt echo as 
a function of time, for various values of the anisotropy 
A and fixed coupling strength e = 0.1. In the ferro- 
magnetic zone outside of the critical region (A > 1) £ 
does not decay at all, as a consequence of the ground 
state full polarization. In the critical region and for long 
times, the Loschmidt echo decay is proportional to the 



modulus of A and for A > —1 it slows down until it is 
completely suppressed in the perfectly anti-ferromagnetic 
regime (A —> —00). The short time decay is again Gaus- 
sian at short times and the rate a is shown in Fig. [12] 
for various sizes of the bath. Wc notice two qualitatively 
different behaviors at the boundaries of the critical re- 
gion: at A = +1 there is a sharp discontinuity, while at 
A = — 1 the curve is continuous. In the critical region 
— 1 < A < 1 the initial decay rate a is constant and 
reaches its maximum value due to the presence of low 
energy modes, while in the ferromagnetic phase A > 1 
it is strictly zero. In the figure, finite size effects are 
evident: indeed, contrary to the AF-model, the decay 
rate a changes with the bath size A^ outside the critical 
region. However, as the system approaches the thermo- 
dynamic limit (A^ 00), the dependence on A^ weakens 
and, while at A = — 1 the curve Q!(A) appears to remain 
continuous, its first derivative with respect to A tends 
to diverge. On the contrary, the ferromagnetic transition 
A = 1 is discontinuous independently of A. 
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FIG. 12: Scaling of the decay rate a at short times for an 
open ended XXZ-Heisenberg chain as a function of A. The 
various symbols are for different sizes of the chain: A'^ = 10 
(black circles), 20 (red squares), 50 (green diamonds), 100 
(blue triangles). The coupling strength between the qubit 
and the chain is kept fixed at e = 0.1. 



C. Decoherence and entanglement 

We propose here to establish a link between decoher- 
ence effects on the system and entanglement inside the 
environment. This study can be justified by the fact that 
decoherence properties of the qubit system seem to be 
intrinsically related to quantum correlations of the bath, 
as the proximity to critical points reveals. On the other 
hand, entanglement quantifies the amount of these cor- 
relations that do not have classical counterparts, and it 
has been widely studied in the recent years, especially 
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in connection with the onset of quantum phase transi- 
tions [3i,[3lj3i,[3i|. 

In Ref. [ill it was shown that two-party entangle- 
ment in the environmental bath of a central spin model 
can suppress decoherence; this effect has been explained 
as a consequence of entanglement sharing, and it was 
supposed to be common to any system whose environ- 
ment maintains appreciable internal entanglement, while 
evolving in time. We now characterize a more complex 
case of system-plus-bath coupling, given by Eq. IH), with 
a richer structure in the ground state entanglement which 
suggests the following picture, valid at short times for 
Ji < 1. 

We expect that when the decay of coherences is 
quadratic, only short range correlations in the bath are 
important, therefore it seems natural to relate the short- 
time decay rate a to the two-site nearest neighbor en- 
tanglement of the chain. All the information needed to 
analyze bipartite entanglement between two spins inside 
the chain, say i and j, is contained in the two-qubit re- 
duced density matrix py, obtained after tracing out all 
the other spins. We use the concurrence C{\i — j\) as 
an entanglement measure for arbitrary mixed state pij of 
two qubits defined as 



C(N - j\) = niax(Ai - A2 - A3 - A4, 0), 



(34) 



where Ai are the square roots of the eigenvalues of the 
matrix PijPij, in decreasing order; the spin flipped matrix 
of Pij is defined as pij — {cry (X) ay) p*j {ay (g) CTy), and the 
complex conjugate is computed in the canonical basis 
{| TT), I Ti), I it), I ii)} • The entanglement of formation 
of the mixed state pij is a monotonic function of the 
concurrence. 

We start our analysis by considering again the XY spin 
bath: the concurrence C'{k) in terms of one-point and 
two-point smn-correlation functions can be analytically 
evaluated [il],!!^]. As long as 7 ^ 0, this system belongs 
to the Ising universality class, for which it has been shown 
that the concurrence between two spins vanishes unless 
the two sites are at most next-to-nearest neighbor [sgI . 
[37j . The nearest- neighbor concurrence C(l) presents a 
maximum that occurs at A = Am > 1, and it is not 
related to the critical properties of the model. For A > 
Am it monotonically decreases to zero, as the ground state 
goes towards the ferromagnetic product state. For A < 
Am the concurrence decreases and precisely reaches zero 
at Xk — -I- 7)(1 — 7), where the ground state exactly 
factorizes |43l |. 

In the proximity of the quantum phase transition, the 
entanglement obeys a scaling behavior |36i] ; in particular 
at the thermodynamic limit the critical point is charac- 
terized by a logarithmic divergence of the first derivative 
of the concurrence: d\C{l) ~ In |A— Ac|. A similar diver- 
gence has been found for the first derivative of the decay 
rate a of the Loschmidt echo (see the inset of Fig. [4]). In 
Fig. [13] we analyzed the dependence of a as a function 
of the nearest-neighbor concurrence C(l). The behavior 
a(C) is not monotonic, since a monotonically decreases 
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FIG. 13: Dependence of the Loschmidt echo decay rate a on 
the nearest-neighbor concurrence C(l) in the XY-model, for 
e = 0.25 and different anisotropy values: 7 = 1 (upper left), 
0.5 (upper right), 0.25 (lower left), 0.1 (lower right). Data 
have been evaluated for an A'^ = 301 spin chain. 



with A (see Fig. [SId), while C(l) has a non- monotonic 
behavior. Indeed, for A > Am and A < Xk the two quan- 
tities are correlated, while they are anti-correlated for 
Am > A > Xk- 



0.01 




FIG. 14: Decay rate q as a function of the nearest-neighbor 
concurrence C(l) in the Heisenberg model for A*' — 100, e = 
0.1. Inset: nearest- neighbor concurrence as a function of the 
anisotropy A. The two dashed lines isolate the critical region. 



A similar scenario has been found in the case of the 
Heisenberg chain as a spin bath, reported in Fig. [141 
As for the Ising model, the behavior of a{C) is not 
monotonic, following the concurrence behavior: Indeed 
C(l) is zero for A > 1 (in the fully polarized ground 



13 



state), it monotonically increases while decreasing A up 
to A = — 1, and then starts to decrease again, vanishing 
in the perfect anti- ferromagnetic hmit A — > — cx) [44]. In 
conclusion, the behavior a = a{C) is monotonically in- 
creasing in the anti-ferromagnetic region, is constant at 
the criticality, and is strictly zero for the ferromagnetic 
phase (see Figs.lTTb. lT^ . resulting in a direct correlation 
of a and C(l) in the anti-ferromagnetic phase. 



RESULTS: THE MULTIPLE-LINKS 
SCENARIO 



In this Section we study the case in which the system S 
is non-locally coupled to more than one spin of the bath 
E, namely we choose to 7^ 1 in the interaction Hamilto- 
nian Eq. In general, given a number to of links, there 
are many different ways to couple the two-level system 
with the bath spins. In the following we will consider two 
geometrically different setups: a spin-symmetric configu- 
ration (type A) where the qubit is linked to some spins of 
the chain which are equally spaced, and a non-symmetric 
configuration (type B) where the spins to which the qubit 
is linked are nearest neighbor. We will also suppose that 
the interaction strength e between the qubit and each 
spin of the bath is kept fixed and equal for all the links. 
A schematic picture of the two configurations is given 
in Fig. 1151 We will present results concerning an en- 
vironment constituted by an Ising spin chain with pe- 
riodic boundary conditions. Notice that this system- 
environment coupling can induce the emergence of criti- 
cality in E even for small cqmDing e <C J in the central 
spin scheme, i.e. m = N |T^. Indeed, in this 

case, a perturbation e results in a change of the external 
magnetic field of the chain from A to A -I- e. If the environ- 
ment is characterized byA<l, ate>l — A the coupling 
drives a quantum phase transition in the bath p^ . 
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Configuration B 
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FIG. 15: DiflFerent setups for the coupling of the two-level 
system (red spin) with the spin chain. The qubit is linked 
to the spins of the chain depicted in grey. Configuration A 
shows a star-symmetric setup, while configuration B is a non- 
symmetric packed setup. As a pictorial example, here we set 
iV = 18, m = 6. 



1. Short-time behavior 

At short times the Loschmidt echo exhibits a Gaussian 
decay, as in Eq. (f22|) . for both configurations. For the 
setup A the decay rate a, far from the critical region, 
scales as 

a oc TO e^, (35) 

as it can be seen from Fig. [T6b . The scaling with the 
number of links to is a consequence of the fact that the 
short-time behavior is dominated by the dynamics of the 
environment spins close to those linked to the qubit (see 
also Sec. IIV C[) . Therefore, if the linked spins are not 
close to each other, they do not interact among them on 
the short time scale Jt <^ 1. Near the critical point this 
picture is not valid, since long-range correlations between 
the spins of the bath become important, even at small 
times. In this region indeed the scaling a '--^ to is less 
appropriate, as it can be seen in the inset of Fig ll6h . 
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FIG. 16: Decay rate a rescaled with respect to the number of 
links m, in a star symmetric configuration (type A) on the left, 
and in a non-symmetric configuration (type B) on the right. 
The environment is an A'^ = 300 spin Ising chain with periodic 
boundary conditions. The interaction strength e is kept fixed 
for all the links between the qubit and the environment spins: 
e = 0.25. Various symbols stand for different values of m: 1 
(plus), 2 (triangles right), 3 (triangles down), 5 (triangles left), 
10 (triangles up), 30 (diamonds), 100 (squares), 150 (circles), 
300 (dashed line, right figure). The inset on the left figure 
shows a magnification of the same plot, centered in proximity 
of the critical point Ac = 1. 



The scaling of Eq.(l35|l does not hold any more in the 
setup B, as in this case collective modes influence the 
dynamics of the system even at short times, since the 
spins linked to the central spin are close to each other 
and they are not independent any more. This is clearly 
visible in Fig.fTBb. We notice also that, as m/N increases, 
a tends to remain constant for A < 1 and then decreases 
for A > 1. In the limiting case to = TV, a is positive 
and strictly constant for A < Ac; the first derivative dxa 
presents a discontinuity at Ac, showing a divergence from 
the ferromagnetic zone A — > 1"*". 
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2. Long-time behavior 



provided that 



In this Section we concentrate on the setup B, since the 
configuration A is less interesting: indeed for long times 
the Loschmidt echo behaves as if the qubit was coupled 
to an environment with a smaller number of spins. In 
Fig. [17] we plotted the asymptotic value of the Loschmidt 
echo Coo for the setup B. We observe that, as m increases, 
the coherence loss enhances, and the valley around the 
critical point Ac deepens and gets broader. Notice also 
that, when approaching the central spin limit, the asymp- 
totic value Coo reaches values very close to zero, even 
far from criticality. This situation is completely differ- 
ent from what occurs in the single-link scenario where, 
away from criticality, the Loschmidt echo remains very 
close to one even for non negligible system-bath coupling 
strengths (e.g., for e = J/4 we found Coo ^ 0.99, as it 
can be seen, for example, from Figs. [21 (H). 



L 



X < 1 and A + e > 1 



(37) 




FIG. 17: Asymptotic value of the Loschmidt echo for the 
setup B; the environment is a 300 spin Ising chain, with e — 
0.25. Various symbols stand for different values of m: 1 (plus), 
2 (triangles right), 5 (triangles down), 10 (triangles left), 30 
(triangles up), 100 (diamonds), 150 (squares), 300 (circles). 



A. Strong coupling to the bath 

Under certain conditions decoherence induced by the 
coupling with a bath manifests universal features: In 
Refs. [ij, [13] it has been shown that when the coupling 
to the system drives a quantum phase transition in the 
environment, the decay of coherences in the system is 
Gaussian in time, with a width independent of the cou- 
pling strength. In particular, for a central spin coupled 
to an iV-spin Ising chain, the Loschmidt echo in Eq. (|C6|) 
is characterized by a Gaussian envelope modulating an 
oscillating term: 



The oscillations are not universal, but the Gaussian 
width Sf^ depends only on the properties of the envi- 
ronment Hamiltonian, in particular it is independent of 
the coupling strength e and of the transverse magnetic 
field A, while it is proportional to the number TV of spins 
in the bath. The case of the central spin model is illus- 
trated in Fig. [T5b . where the different curves stand for 
various values of e. 




C{t) = |cos(et)|^/2 



(36) 



FIG. 18: Loschmidt echo for the Ising model with N = 300 
spins as environment, A = 0.5. The various panels are for 
different numbers of spins m coupled to the qubit system and 
for various geometries: a) m = 300 (central spin system); b) 
m = 3, configuration B; c) m = 30, configuration A; d) m = 
30, configuration B. Different curves correspond to various 
values of the coupling strength: e = 2 (full black), 5 (full 
red), 10 (full green), 20 (dashed blue), 40 (dotted magenta). 



This universal behavior for strong couplings is present 
in more general qubit-environment coupling, different 
from the central spin limit, provided that the same con- 
ditions in Eq. ([57)) on A and e hold. In the Ising bath 
model, we checked that it remains valid also for a num- 
ber of system-bath links m ^ N and for different ge- 
ometries, as it can be seen in Fig. [T8| (panels b-d). The 
fast oscillations remain proportional to e, since they are 
a consequence of the interaction oc ^TzCJ^ between the 
qubit and the spins of the bath. Instead in general the 
envelope is Gaussian only for small times, but its width 
is independent of e. 

We found that, for the non-symmetric configuration B, 
S]^ depends only on the number of links m, being pro- 
portional to it. No dependence from the perturbation 
strength e, the size A'' of the bath, or the transverse mag- 
netic field A has been observed, provided that m > 10, 
as it can be seen from Fig. [T^ . Indeed in the setup B 
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the qubit behaves hke in the central spin model: in the 
regime of strong coupling, the effect of the environment 
spins not directly linked to the qubit can be neglected 
since they interact weakly, as compared to the system- 
bath interactions. Therefore a system strongly coupled 
to m neighbour spins of an A^-spin chain decoheres in 
the same way as if it was coupled to all the spins of an 
m-spin chain. 




m m 



FIG. 19: Width of the Gaussian envelope S% as a function 
of the number of links m between the system and the bath, 
for different values of the transverse magnetic field: A = 
(circles), 0.1 (squares), 0.2 (diamonds), 0.5 (triangles up), 0.7 
(triangles down), 0.9 (triangles left). The bath is modeled as 
an A'^ = 200 Ising chain with periodic boundary conditions; 
the coupling strength is kept fixed and equal to e = 80. The 
two panels correspond to: a) star-symmetric configuration 
(type A); b) non-symmetric configuration (type B). Dashed 
black curves represent a linear behavior S% ~ m and are 
plotted as guidelines. 



The decoherence effects are richer in the geometry A. 
In this setup the effect of the spins not directly coupled 
to the qubit cannot be ignored, indeed we found that 
depends on the internal dynamics of the bath. In 
particular we found that 

5^ - m A^ (38) 

where A is the transverse magnetic field (see Fig. [TOk and 
Fig. [20)1 . Notice that, as it can be seen from Fig. [T9k . 
the case A = is rather peculiar, since the system does 
not decohere at all, unless it is uniformly coupled to all 
the spins of the environment. This is due to the fact 
that an Ising chain in the absence of a transverse field 
is not capable of transmitting quantum information: in 
this system the density matrix of two non-neighbouring 
spins evolves independently of the other spins. 



VI. CONCLUSIONS 

In this work we analyzed a model of environment con- 
stituted by an interacting one-dimensional quantum spin- 
1/2 chain. We have shown that this type of baths can 



10^ FT 




FIG. 20: S% as a function of the transverse field A in an 
= 200 spin Ising bath, for different numbers of system- 
bath links in the configuration A: m = 2 (circles), 4 (squares), 
10 (diamonds), 25 (triangles up), 50 (triangles down), 100 
(triangles left), 200 (triangles right). The coupling strength 
is kept fixed and equal to e = 80. The case m = 200 (triangles 
right) corresponds to the central spin model. 



be experimentally engineered in a fully controllable and 
tunable way by means of optical lattices, and we pro- 
vided a detailed description of the pulse sequence needed 
to realize a pure dephasing model for a qubit system cou- 
pled to a spin bath. The open quantum system model 
we proposed here is simpler to implement in an experi- 
mental setup, with respect to the widely studied central 
spin model in which the system is uniformly coupled to 
all the spins of the environment. Nonetheless, we have 
analyzed the crossover from the single link to the central 
spin model by taking in account different geometries and 
different number of links, recovering, in the central spin 
limit, known results present in literature, and character- 
izing new behaviors in between. 

The coherence loss in the system has been quantified 
by means of the Loschmidt echo; the rich structure of 
the environment baths allowed us to analyze the impact 
of the different phases and of the distance from critical- 
ity, and more generally to study the connection between 
entanglement inside the bath and decoherence effects in 
the system. 

First we considered the case of one link between the 
qubit and the chain. We showed that at short times 
the Loschmidt echo decays as a Gaussian, while for long 
times it approaches an asymptotic value which strongly 
depends on the transverse field strength. Criticality 
emerges in both short- and long-time behaviors: the XY 
spin bath is characterized by the first derivative of the 
Gaussian width with respect to the field strength that 
diverges logarithmically at Ac, while the saturation val- 
ues of the Loschmidt echo at long times exhibits a sharp 
cusp in correspondence to the critical point. The de- 
cay rate at short times due to a Heisenberg spin bath 
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is maximum and constant throughout the whole criti- 
cal region. We also showed that, contrary to the central 
spin model, if the environment is not critical, the decay 
of coherences in the system is independent of the bath 
size. Subsequently we related decoherence properties of 
the system with quantum correlations inside the environ- 
ment: we studied the Gaussian decay rate at short times 
as a function of the concurrence of two nearest neighbor 
spins in the bath. We then analyzed a non-local system- 
bath coupling in which the qubit is connected to the chain 
with multiple links, by considering two different setups: 
a spin-symmetric and a non-symmetric one. Finally we 
reported the existence of a coupling-independent regime, 
with the only hypothesis of having a strong system-bath 
coupling. 
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APPENDIX A: DMRG APPROACH 

An analytic solution for a spin coupled to a Heisen- 
berg XXZ-model (i.e. {A 7^ 0, A = 0} in Eq. ^) is 
not available, therefore a numerical approach is required. 
We resort to the recently developed time-dependent Den- 
sity Matrix Renormalization Gro up ( t-DMRG) with open 
boundary conditions [H, H^, H^, IZJ), in order to calcu- 
late the overlap between the ground state |G) of Hg and 
the time evolution of the same state under the Hamilto- 
nian He (see Eq. ©). 

First, the static DMRG algorithm allows us to eval- 
uate the ground state \G) of Hg. This algorithm in its 
first formulation 45| is an iterative numerical technique 
to find the ground state of a one-dimensional system con- 
stituted by sites which possess local and nearest-neighbor 
couplings, therefore it is well suited for systems like the 
Hamiltonian in Eq. ^ . The key strategy of the DMRG 
is to construct a portion of the system, called the system 
block, and then recursively enlarge it, until the desired 
size is reached. At every step the basis of the corre- 
sponding Hamiltonian is truncated, so that the size of 
the Hilbert space remains manageable as the physical 
system grows. The truncation of the Hilbert space is 
performed by retaining the eigenstates corresponding to 



the D highest eigenvalues of the block's reduced density 
matrix. 

The t-DMRG is then subsequently used, in order to 
simulate the dynamics of \G) under the Hamiltonian Tig. 
The t-DMRG algorithm |46!] is an extension of the static 
DMRG, which follows the dynamics of a certain state 
lip) under the Hamiltonian 7i of a nearest-neighbor one- 
dimensional system. The time evolution can be imple- 
mented by using a second-order Suzuki- Trotter decom- 



position [4^ of the time evolution operator U 



i' 



(Al) 



Here we have first written the Hamiltonian as Tip = T+Q 
where !F contains the on-site terms and the even bonds 
while Q is formed by the odd bonds; n — t/ 5t is the time 
expressed in number of trotter steps. We also keep track 
of the state \G) in the new truncated basis of 6"'^^=* |G), 
so that we can straightforwardly evaluate the overlap in 
Eq. ([9]) . During the evolution the wave function is chang- 
ing, therefore the truncated basis chosen to represent the 
initial state has to be updated by repeating the DMRG 
renormalization procedure using the instantaneous state 
as the target state for the reduced density matrix. 

The results presented in this paper concerning Heisen- 
berg spin baths, have been obtained by performing time 
evolution within a second order Trotter expansion, with a 
Trotter slicing J St ^ 10~^ and a truncated Hilbert space 
of dimension D = 100. The evaluation of the Loschmidt 
echo with the t-DMRG is much more time- and memory- 
expensive than the analytical approach of Eq. pT|) . lim- 
iting the bath size in our simulations up to ~ 10^ 
spins. 



APPENDIX B: FERMION CORRELATION 
FUNCTIONS 

We provide an explicit expression for the two-point 
correlation matrix {r)ij = (\E'|\E'j) of the operators = 



(c| . . . c]y ci . . . Cat) on the ground state of the system 
Hamiltonian Hg. The 2N x 2N matrix r is written in 
terms of the Jordan Wigner fermions {cfc, cj,} as: 

■ (Bl) 



N 



Therefore it is sufficient to express {c^, c^} in terms of 
the normal mode operators {Vk^K vlf^^ } which diagonal- 
ize the Hamiltonian Hg, since < rjj^^ijlf^'^ >= djk, the 
other expectation values of the r/'s on the ground state 
being zero. By inverting Eq. (|18p we get 

c = g'^ • 7? + h"^ • , (B2) 

from which it directly follows that 



h(g)T h(s) h(s)T g(g) 
g(g)T h(g) g(g)T g(g) 



(B3) 
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The superscript stands for the change-of basis- 
matrices g and h relative to the Hamiltonian Hg. 

The last ingredient for the evaluation of the Loschmidt 
echo, Eq. (PT|) . is the exponential e'*-'*. We introduce 
the vector pt = (t^J^H . . . t^^H ^(s) . . . , so that * = 
IJt r, where 



U = 



We can therefore rewrite Eq. l|20|) as 



(B4) 



K = ^*^c* = ir^ucu^r = ir^Dr, (B5) 

where D is a 2N x 2N diagonal matrix, whose elements 
are the energy eigenvalues of Tig and their opposites: 



D 



-E^'') 



(B6) 



It then follows that C = U'l' D U, from which one can 
easily calculate the exponential e**^*. 



The corresponding single quasi-excitation energy is given 
by4^)^ffe(e): 



/27rfc\ , 



2 ■ 2 /27rfcN 
7^sm (— j 



(C3) 



The Hamiltonian Tig can be diagonalized in a simi- 
lar way; the corresponding normal-mode operators are 
connected to the previous ones by a Bogoliubov transfor- 
mation: 



?41 = cos(afe) ry^^^ - i sin(afe) r?J]t ^ (C4) 
with ak = MO) - 9k{e)]/2. 

The spin chain is initially in the ground state \G) of 
Hg; this state can be rewritten as a BCS-like state: 



|G) = n [cos(a,.) -^sin(afe)r;i'=H^(_-)t] (C5) 

/£>0 



APPENDIX C: CENTRAL SPIN MODEL 

If the qubit is uniformly coupled to all the spins of 
the chain, the effect of the interaction in Eq. (|4]) is sim- 
ply of renormalizing the transverse field strength in the 
bath Hamiltonian: for A = both Hamiltonians Tig and 
Tie correspond to an Xl"-model with anisotropy 7, uni- 
form couplings J and local magnetic field A and A -I- e 
respectively. They can be diagonalized via a standard 
JWT, followed by a Bogoliubov rotation [3l|. The nor- 
mal mode operators that diagonalize Tie, satisfying the 
fcrmion anti-commutation rules, are given by 



(e) e 



-2iTijk/N 



j 



N 



(e) + 



(CI) 



where the coefhcients u 
depend on the angle 

^fc(e) = arctan 



cos 



{ek/2), 4"^^ = sin(0fc/2) 



-sm{2TTk/N) 



cos(27rfc/7V) - (A -He) 



(C2) 



where \Ge) is the ground state of Tie- This expression 
allows one to rewrite the Loschmidt echo C{t) in Eq. ^ 
in a simple factorized form: 



N/2 



C^=N{t) = n [1 - sin2(2afc) sin2(ffei)] 



(C6) 



fc=i 



Eq. (jC6[) provides a straightforward formula in order 
to calculate the Loschmidt echo for a central spin coupled 
uniformly to all the spins of the bath. Nonetheless this 
model, although in some circumstances it can reveal the 
emergence of a quantum phase transition in the environ- 
ment [H, [13, [3 , 'ioss '^ot provide an effective physical 
description of a standard reservoir, since the coupling 
with the system is highly non-local and it can drive the 
evolution of the bath itself. However it is useful through- 
out the paper, to test the convergence of our results to 
the ones presented in [ll, [13, [l^ , in the limit m ^ N. 
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